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Abstract 

We associate a fractal in RP" to eacli vector basis of R"+^ and we 
study its measure and asymptotic properties. Then we discuss and study 
numerically in detail the cases n — 1,2,3, evaluating in particular their 
Hausdorff dimension. 



1 Introduction 

In this paper we study an algorithm that takes a basis of and builds, out 

of it, a fractal in RP". We do this by using the following two basic facts: 1. In 
RP", every n + 1 points {pi} which do not lie on the same afhne hyperplane 
determine a partition of RP" in the 2" projective n-simplices having the points 
Pi as vertices; 2. Given a vector basis £ — {ei} in R"+^, we can build n+1 new 
bases £i by fixing the i-th vector and summing it to the n remaining ones. 
Now, consider a vector basis £ of R"+^. Its vectors et projects into n + 1 point|l| 
[ci] € RP" and therefore determine a partition of RP" in projective n-simplices 
as in point (1). We denote by S{£) the one the point [ei + 62 + 63]. Next, 
consider the n+1 projective n-simplices S{£i) corresponding to the bases £i 
defined in (2). Their union U^+/S'(f,) can be thought as the difference between 
S{£) and the interior of the projective polytope Z{£) (the body of £) having 
the points [e^ + ej] as vertices. By repeating recursively this step on the S{£i), 
we end up building a (n + l)-ary tree of bases T{£) = {£1}, where I = ii . . .ik 
is a multiindex, and fractal F{£) whose points are the ones left inside S{£) after 
removing all the bodies Z{£j). 

Topologically, this fractal coincide with the multi-dimensional generalization 
of the Sierpinski triangle [Siel5| , namely the fractal generated by removing from 
a n-simplex S the polytope Z having as vertives the middle points of the edges 
of S. Geometrically though they are different because the vertices [ei + ej] of 
the body Z{£) are closer to the vertices oi S{£) corresponding to the vectors of 
higher Euclidean norm, and even if we start with a basis where all vectors have 
the same norm they will not be anymore so after the first step. 

-'Here and throughout the paper we denote by [e] € RP" the direction of the vector 
e e M"+i. 
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We were motivated to study this fractal by the following two reasons. First, 
this fractal is invariant with respect to a set of n + 1 transformations (more 
specifically, projective diffeomorphisms) but they fail to be an Iterated Func- 
tion System (IFS) because they are not contractions. In particular the ma- 
chinery developed for IFSs cannot be applied to this case and no analytical 
bounds are known for the Hausdorff dimension of F{£) (numerical estimates 
for the cases n = 2,3 can be found in Section [4]). Second, the construction 
leading from £ to F{£) is the n-dimensional generalization of the structure dis- 
covered by the author and LA. Dynnikov in a fractal, in MP^, which describes 
the asymptotic behaviour of the plane sections of the triply-periodic cubic poly- 
hedron C = {4,6|4} [DDO9JI. Indeed one of the results of Novikov's theory of 
plane sections of triply-periodic surfaces [Nov82i |Dyn99| is the following. Con- 
sider a connected triply-periodic surface S which divides in two components 
which are equal modulo translations. Then there exists no bundle of paral- 
lel planes whose intersections with S are all compact. The asymptotics of the 
open (i.e. non-compact) sections, as function of the direction of the bundle 
of planes, are described by a "labeled cut-out fractal" in the projective plane 
(see |De 03[ IDe 06] for a few other concrete examples). Here by labeled cut- 
out fractal we mean a fractal F which is obtained by removing, from an initial 
region, a sequence of closed sets {Zi} whose interiors is pairwise disjoint and 
such that to each Zi it is associated an element bi of some set of "labels" B. 
In the particular case of Novikov's theory above, the set of labels is the set of 
all indivisible triples of integers and the geometrical meaning of the fractal is 
the following. Let d be any vector in M'^ whose direction [d] belongs to some 
Zi labeled by bi. Then the open sections arising by cutting the surface S with 
planes perpendicular to d are strongly asymptotic to a straight line whose di- 
rection is given by the "vector product" d x b. Going back to the surface C, it 
turns out that its corresponding fractal Fc essentially coincides with the frac- 
tal F{£) c MP^ associated to the basis £ = {(1, 0, 1), (0, 1, 1), (1, 1, 0)} of K^. 
We conclude this long digression by pointing out that the Novikov's theory of 
plane sections of triply-periodic surfaces is the mathematical model for the phe- 
nomenon of the anisotropic behaviour of magnetoresistance in normal metals 
at low temperature and under a strong magnetic field (see |LP60| and |NM03| 
for more details about the physics and the dynamics of this phenomenon). In 
particular Fc encodes the information on the conduction of the electric current 
in a metal having C as Fermi Surface. 

The paper is organized as follows. In Section [2] we define the basic objects 
and prove a few elementary facts about them. In Section [3] we prove that all 
fractals F{£) have zero volume with respect to some natural measure on the 
projective space and study their asymptotic properties, showing in particular 
that they are related with the n-bonacci sequences. Finally, in Section [4l we 
discuss in detail the cases n = 1,2,3 and present numerical results indicating 
that, unlike the Sierpinski case, the Hausdorff measure of F{£) may be non- 



^We discovered later that this fractal had been already considered in the past by 
G. Levitt | Lev93| while studying dynamical systems on the circle. 
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integer even for n = 3. 

2 Structure of the fractal 

Let £ = {ei, . . . , e„+i} bo a vector basis of M"+-'^ and let us call volume of £ 
the Euclidean volume of the (n + l)-simplex of M"+^ naturally associated to it. 
Every of the n + 1 sets £j = {e\, . . . , e^j+i}, j = 1, • • • , n + 1, defined by 



is also a vector basis of R"+^ and has the same volume as £. Repeating re- 
cursively this procedure, we get an inifnite (n + l)-ary ordered rooted tree 
T{£) = {£l^}k(^^ of bases of all with the same volume, with £ as root. 

The multi-index Ik = «i«2 • • • ik-iik describes the steps needed to build the basis 
from the root, namely £i^^ii^ = (f/t.Jifc. 

This tree structure corresponds to the limit process for building a fractal on 
RP". Indeed let £ = {ei, . . . , e„+i} and denote by [e^] G MP" the direction of 
the vector e^. To £ it is naturally associated a projective n-simplex S{£) defined 
in the following way. The n+1 points [ci] G MP" are the vertices of 2" projective 
n-simplices whose interiors are pairwise disjoint and whose union gives the whole 
R-P"; the point p = [ei -I- • • • -I- e„+i] is not a boundary point for any of them and 
we denote by S{£) the one which contains p. Now consider the bases £i at the 
first recursion level of T{£). The vertices of the projective n-simplex S{£i) are 
the points {[ei -|- e^], . . . , [ei_i -I- ei], [e^], [ci+i -|- e^], . . . , [e„+i +ei]}, i.e. S{£i) is 
contained inside S{£), shares with it the vertex [e^] and (part of) all the edges 
coming out from that point and has in common exactly one vertex with each 
other S{£j), j ^ i. 

Let now Fi — iJl^iS{£i). The difference between S{£) and Fi is the interior 
of the projective polytope Z{£) having the n{n + l)/2 points [e-i + Cj], i ^ j, 
as vertices. We call Z{£) the body of £. More generally, let Fk = U\i\=kS{£i) 
be the fc-th level of recursion of the fractal, with Fa = S{£). The set Fk is 
obtained from Fk-i by erasing the interiors of the {n+l)''~^ bodies Zj = Z{£j), 
\J\ = k — 1. The fractal F{£) is then obtained as the hmit F{£) = flfegN^fe- 

Note that we can always find an affine n-plane (i.e. a canonical chart for 
the projective space) inside RP" which contains the entire S{£) and, therefore, 
the whole F{£). From now on then we will consider often S{£) and all the 
S{£i) as an n-simplex inside M". This allows to provide another geometric 
characterization of the algorithm generating the fractal. Indeed the fc-skeleton 
of S{£) is the set of the convex hulls associated to the ("^^) different subsets of 
k elements of £, namely the convex hulls of the sets {[ej^], . . . , [e,^]} where no 
two indices are equal. 

Definition 1. The vector h{£) = Yll=i ^ M"+^ is called the barycenter of 
the n- simplex S{£), where £ = {ei}i=i^... ^n+i- Analogously, the barycenter of its 
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k-face of vertices {[e^J, . . . , [e^^]} is the vector bj{£) = J2'j=i ^-ij '= where 
I = ii.. .ik- By abuse of notation we sometimes call barycenter its direction 
[b] e KP" . It will be clear from the context which one we are referring to. 

Lemma 1. Let £ be a basis o/]R"+^ and fi^ the k-subsimplex of S{£) corre- 
sponding to {e^^, . . . , e^^} C £. Then the projection of the barycenter \bi,,] of 
fi^ from the vertex [cj^.] on the {k — l)-face of fi^. opposite to it, namely the 
one corresponding to the k — 1 vectors {e,!, . . . , ej^} \ {ej^.}, coincides with the 
barycenter of that face. 

Proof. This relation is clearly recursive and therefore it is enough to prove the 
theorem in the case of the n-simplex S{£) and any of its faces. Let us consider 
what happens for the vertex [ei]: the face fi of S{£) opposite to it corresponds, 
in R"+^, to the n-planc fi spanned by the n vectors 5*^^^ = {ej}j^i and the line 
li joining [ei] to [b] corresponds to the 2-plane li spanned by ei and b. Since 
b = J22=i ^k, clearly the only linear combinations belonging to both /i and /i 
are the span of the vector Y^^=2 = b — ei. In other words, the intersection 
between /i and Zi is [^2 + • • • + e„+i], which is indeed the barycenter of £'^^^ and 
similarly for the other vortices. □ 

Proposition 1. The vertices of the body Z{£) corresponding to a basis £ can 
be obtained in the following way: project the barycenter of £ from its vertices to 
its faces and repeat recursively this procedure until the edges are reached. The 
n(n + l)/2 points obtained are the vertices of Z{£). 

Proof. The recursive procedure makes sense because, thanks to the previous 
lemma, we know that the projection of the barycenter on a face via the vertex 
opposite to it coincides with the barycenter of the face. When we reach the 
edges, therefore, we are left with their barycenters, which are clearly the n{n + 
l)/2 points [cj + Cj], i ^ j. □ 

Finally, we provide a third way to describe this fractal. Recall that the 

Sierpinski gasket and its natural multi-dimensional generalization can be seen 
the invariant set of a Iterated Functions Systems (IFS). Similarly, we prove below 
that the fractal F{£) is the invariant set of n + 1 projective diffeomorphisms 
{"fpi}. They do not form however, strictly speaking, a IRS because they are not 
contractions; in particular the Jacobian of each of them is the identity in the 
omonimous vertex [e,] of S{£). 

Proposition 2. The fractal F{£) is invariant with respect to the (n+ 1) projec- 
tive automorphisms ipi of RP" induced by the linear transformations -tjji defined 

by -i/Jiicj) = Ci + Cj, j ^ i, ipiie^) = e^. 

Proof. This is simply a consequence of the fact that the tree itself T{£) is clearly 
invariant under the action of the -tjji, so that the ipt map the set of bodies Z{£j) 
into itself and therefore leave the fractal invariant. □ 

Remark 1. Every body Zi{£), \I\ = k, is the image of the root body Z{£) via 
the map tpi := ipi^ o ■ ■ ■ o tpi^^ . 
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3 Measure and Asymptotics of the fractal 

We start by proving that the volume of F(£) is zero with respect to any measure 
induced on MP" by the Lebesgue measure on any affine n-plane or, equiva- 
lently, with respect to the measure induced by the canonical one on the sphere. 
We begin with a technical lemma: 

Lemma 2. The maximum mk,n of the functions 

n 

fk,n{vi, ■ ■ ■ ,Vn) = ^ n n n 



i=2 i=l i=2 



where fc,7i G N, n > 2, on the n-simplex S with vertices in pq — (0, • • • ,0), 
pi = (1, 0, • • • ,0), • • • , p„ = (0, • • • , 0, 1), is given by 

on 

mk,n = /a;,„(1,0, • • • ,0) = 



(2 + fc)(3 + fc)" 

except for the case n = 2, k = 0, where mo,2 = /o,2(0, 0, • • • ,0) — 1/4 

Proof. A direct computation shows that the derivative of fk,n with respect to 
any w^, i > 1, is negative inside S and therefore the maximum is attained at the 
smallest values for those variables. Then we are left with the function 

f (l + ^i)" 

''"^"'^^ {l + k + v^){2 + k + v^y- 

whose derivative 

(1 + + (n - 2 - + vi) + k[n{2 + wi) - 1 - vi]} 



{l + k + vi)^{2 + k + vi) 



n+l 



in the domain vi g (0, 1) is always positive for fc > while for fc = is always 
positive for n > 2 and always negative for rt = 2. 

Hence fk^n will reach its maximum in the origin when fc = 0, n = 2 and in 
the point (1, 0, • • • , 0) in all other cases. □ 

The following proof is a generalization to any n of the proof provided in [DD09j 
for the case n = 2. 

Theorem 1. The fractal set F{£) is a null set for fi. 

Proof. It is enough to prove the therem for a particular choice of £ = {ei = 
(1,0,-- - ,0,1),-- - ,e„ = (O,-- - ,0,l,l),e„+i = (0, • • • ,0,1)}. With this £, the 
n-simplex S{£) is contained in the affine plane tt = {/i„+i =/= 0} with respect 
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to the homogeneous coordinates [hi : ■ ■ ■ : hn+i]- On tt wc use the canonical 
coordinates {vi,. . . , u„) defined by Vi = hi/hn+i and the measure 

dvi ■ • • dvn 



Each of the ipi maps the whole fractal F = F{£) c <S'(^) into disjoint sets 
^i(F) = F{Si) C so that fiiF) = J2"^d KMF))- With our choice of £ and 
fi both S{£) and the measure (at least close to S{£)) are invariant with respect 
to every permutation of the first n basis vectors and therefore fi{tjji{F)) = 
ij,{tpi{F)),yi = , n, so that fi{F) = n/i(V'i(-F))+/i(V'„+i(F)). By repeating 
this procedure on ipn+i{F) we find that 

SO that jl{F) = n/i(V'i + nfl{ijjn+i o ipi (F)) + (l{i>n+i o i>n+i {F)) and finally, 
by recursion, 

oo 
fe=0 

since limfe^oo = 0. 

Wc will now show that fl{ipn+i ° < with J2T=o'^k^^ < -'■/'^' 

which leads immediately = 0. Note indeed that, with this particular choice 
of the basis, the action of the tpi on the corresponding homogeneous coordinates 
is given by 

'il)i{[hi : h2 ■■■■■■■ hn ■■ hn+i]) = : h2 : ■ ■ ■ : hn : 2hn+i - hi] 



ipn{[hi : h2 ■■■■■■■ hn : hn+i]) = [hi : h2 ■■ ■ ■ ■ ■■ hn+i : 2ft„+i - hn] 
^'ipn+i{[hi : h2 ■■■■■■■ hn ■■ hn+i]) = [hi : h2 ■■■■■■■ hn ■■ J2"^i hi] 

Since the fractal is invariant with respect to the projective transformation 

n+l 

R{[hi : h2 h„: /i„+i]) = [2/i„+i - ^hi : h2 hn : hn+i] 

i=l 

corresponding to the exchange of the vectors ei and e„+i, we can replace f/'n+i ° 
iPi{F) with V^+i o V-i o R{F). 
Then 

n+l 



ipi{R{[hi : h2 : ■ ■ ■ : hn : hn+i])) = [/i„+i : h2 hn : ^ hi] 

and finally 

V'n+i (V'l {R([hi : h2 : ■ ■ ■ : hn ■■ j = [hn+i : h2 hn : ^ hi+k ^ hi] . 



i=l 



n+l n+l 



1=1 



A generalized Sierpinski gasket in 



In the chart Vi = hi/hn+i the map fk — '4'n+i o "01 ° is represented by 

TL n 

fk{vi,--- ,Vn) = {1/D,V2/D,--- ,Vn/D), D ^ 1 + k{l . 



A direct computation shows that the Jacobian of is given by 

1 



SO that Al^f'-*) < for 



max 



As shown in Lemma [2] 



max 



2" 



(2 + fc)(3 + fc)'' 



with the sole exception of the case n = 2, fc = 0, in which case Cg^-* = 1/4. If 
n = 2 then, as aheady shown in [DD09] , 

2^ 

^ (2 + fc)(3 + fc)2 36 3 



fc=0 fe 

In the n > 2 case instead we use the fact that 
„(») ^ 

k=0 



- -TT-^ ~ 0.45 < - . 



fe=0 
2n-l 

3" 

2n-l 



(2 + fc)(3 + /c) 



on— 1 

< + 2" 



n Qn 



da; 



(2 + a;)(3 + a;)" 



2" 



fc=i 



dx 



(3 + x)'= 



By Taylor's expansion theorem applied to log(l — a:) we know that there exist 
a ^ e (0, 1/3) such that 



k 

so that finally 



1 



< 



n ^ 



1 



2 ^ fc3'= (n+l)3"+i(l-£)"+^ ^ fc3'= (n + l)2 

A.;— 1 / \ ^/ k—l 



V,("),2"i 1 1 

^ 3"4^n^+2(n+l) 
fc=o ^ ' 



n+l 



IS 



It is easy to verify that the analytical function g{x) = {'^)^{^ + 2(x+i) 

bigger than h{x) = 1/x for x > 4, which proves that J2T=o'^k^^ ^ V" for ^-ll 
n > 4. 
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We complete the proof by verifying the case 71 = 3 by a direct computation: 

00 00 „3 . 

fro fr(2 + fc)(3 + fc)3 3 ''^^ 3 

□ 

Next Corollary will be used later to justify one of the numerical methods we 
used to evaluate the box-counting dimension of the fractal. An illustration of it 
can be found in Fig. [D 

Corollary 1. The fractal set F{£) is contained in the set of accumulation points 
of the set of barycenters. In particular, the closure of the set of the barycenters 
is equal to the union of F{£) with the boundaries of the bodies Z{£j), £j e T{£). 

Proof. Since F(£) has zero measure it cannot contain any open set. In other 
words, every open set inside S{£) either is contained inside a body Z{£j) for 
some multi- index / or contains one of them. Let p G F{£). Then any open 
neighborhood of p is not contained inside a body and therefore contains one. 
Inside every body lies a barycenter and so F{£) is contained in the closure of 
the (countable) set of barycenters. □ 

In order to study the asymptotics of the fractal it is convenient to pose the 
following definition: 

Definition 2. We call "section" of an infinite tree T a sequence {^ijigN C T 
such that each element t„ (except for the first) is child of its antecedent 

In the Sierpinski case the asymptotics properties do not depend on the par- 
ticular section but in case of F{£) they do. E.g. consider an edge t of T{£), i.e. 
a section t = {t^} defined by t^ = £jj . . . j for some index i. Then the volume 

k 

of the simplices S'(tfc) decreases polynomially with fc, while in the Sierpinski 
case they always decrease exponentially. Below we study the sections where the 
volumes grow faster. They are related to n-bonacci sequence, namely sequences 
whose fc-th element is equal to the sum of the previous n ones, thanks to the 
next proposition: 

Proposition 3. Let A be the set of bodies having non-empty intersection with 
Zi{£) and B the set of indices of the hyperfaces of S{£) (labeled after the index of 
the vertex opposite to it) having points in common with Zi(£). Then b{Zj{£)) = 
Y^zeA K^) + J2keB ^k, where bk = ^^ - nek- 

Proof. Since the fractal is invariant under the ^pi and they are induced by linear 
transformations, it is enough to prove this property for the barycenter b of the 
root cut-out polytope Z = Z{£), which cuts all faces of = S{£). Using the 
^p~^ it is easy to determine that the body corresponding to the k-th face of S 
is the body of the simplex of vertices 



{[ei - Ck], . . . , [ck-i - Sk], [efe], [efc+i - e^], . . . , [e„+i - e^]} 
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Figure 1: Detail, in the — 1 chart, of the first few bodies correspond- 
ing to a Tribonacci section starting from the root of the tree T{£), where £ = 
{(1, 0, 1), (0, 1, 1), (0, 0, 1)}. The first barycenter bi — (1,1,3) of the section is not 
shown. The next five ones, whose projection on RP^ is shown above, are &12 — (1, 3, 5), 
6123 = (3,5,9), 61231 = (5,9,17), 612312 = (9,17,31) and 6123123 = (17,31,57). The 
centers of the bodies of the section lie on a smooth "Tribonacci projective spiral" 
drawn above which is winding about (l/as, l/ai) — (.296, .544). 



and therefore its barycenter is the vector bk — Y!h=i ~ '^^fe- Now it is easy to 
verify that 

n+1 n+1 /n+1 \ n+1 n+1 n+1 

^ = ^ ^ ei - nefe = (n + 1) ^ - n ^ = ^ 6 . 

fc=l fe=l \i=l / i=l k=l i=l 

□ 

This result suggests the following interesting way of building sections of a 
tree T{£) for a basis £ of M"+^. Pick any element ti — and continue the 
section recursively by taking = fii,ii±i,- - ,ii±j, where all indices are meant 
modulo {n + 1). By construction, the body of i^_|_2 touches the bodies of all of 
the previouselements of tf and therefore its barycenter is given exactly by the 
sum of the barycenters of their bodies, and the same happens for all remaining 
terms tf , i > n + 2. We call Fibonacci sections this particular kind of sections 
because the sequence of the corresponding barycenters is a n-bonacci sequence. 
Fibonacci sections of T{£) are relevant for two reasons: 1. they represent the 
sections with faster growth hich barycenters norms grow faster; 2. they provide 
a way to get explicit expressions for points in F{£). 
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Theorem 2. Barycenters of 's bodies in a Fibonacci sequence grow in norm 
as a'' , where a is the {n+ \ )-honacci number (i.e. the highest module root of the 
equation a;"+^ = a;" + • • • + 2; + Ij. This is the highest growth rate for barycenters ' 
norms on a section ofT{£). 

Proof. We can assume without loss of generality that Ci is the canonical ba- 
sis for since asymptotics will not change under the action of a single 
invertible linear transformation, and we can prove the result using the norm 
||t;||i = Y^^=i because in finite dimension all norms are equivalent. 

It is well known that the fc-th term, fc>n-|-l,ofa(n-|- l)-bonacci sequence 
can be expressed as a linear combination with constant coefficients of the k- 
th powers of the n + 1 complex roots of the {n + l)-bonacci equation = 
x^~^ + • • • + X + 1. The highest module root is known to be real and it is 
called (n + l)-bonacci constant. Asymptotically only the highest module root 
is relevant and this proves the first part of the theorem. 

Now, assume that up to the n-th recursive step it happens that at each 
step k the bodies with higher baricentric norm are the ones built starting from 
Z{£) and belonging to a Fibonacci section: then at the following recursive step 
the bodies with higher barycentric norm are exactly the ones which continue 
those Fibonacci sections. Indeed no body can touch more than one body from 
each tree level since bodies corresponding to the same level belong to distinct 
simplices; hence at the {k + l)-th level the bodies' barycenters of the mem- 
bers of those Fibonacci sections are obtained by summing of the highest norm 
barycenters and the components are all positive, so their norm is the biggest 
achievable. □ 

Remark 2. Proposition]^ applied to Fibonacci sections, grants that the limit 
point of a Fibonacci section must belong to F{£). Consider for example the 
Fibonacci sequence generated by 

6_„ = (1, • • • , 1, 1 - n), • • • , feo = (1 - 1, • • • , 1), = (1, • • • , 1) 

In this case all components follow the very same sequence but the component j is 
shifted by one with respect to the component j + 1 for j = 1 , • • • , n, namely 5^ = 
^fcl^i- The last component fe^"*"^ has "initial conditions" b"il^^ — —n, fe"^!*^! = 1, 
• • • , b^^^ — \, so that the first terms of the sequence are ~ 1, 62^"'^ = 

n -|- 1, b^'^^ = 2n + 1 and so on. Since the k-th term of a n-bonacci sequence 
behaves asymptotically like , in RP" the sequence of the corresponding points 
converges to {\ : a : ■ ■ ■ : a"). 

The following theorems shows that edges and Fibonacci sections are respec- 
tively the slower and faster sections with respect to volumes' growth. 

Theorem 3. Let £ be a basis o/R"+^, T{£) ~ {ti} its tree of bases and B{£) = 
{5/} the corresponding tree of barycenters. Then there exist real constants A, B 
such that 

A\I\ < \\bi\\ < Ba\^\ 

for all multiindices I . 
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Proof. We can prove without loss of generality the theorem by fixing the basis 
as the canonical basis of R"+^ and the norm as the maximum norm ||w||oo = 
max \vi\. 

As shown in Theorem [21 the biggest barycenters at every level k are those 
belonging to a Fibonacci section starting by the root element of the tree; the ex- 
plicit expression for those sections, modulo permutations, is bk = {ak-n-i, ■ ■ ■ , o-k), 
k > n + 1, where au — J2"=i ^i'^i^ the ai are the root of the {n + l)-bonacci 
equation and Ai = \/Hj^i(ai — aj). We order the roots so that ai = a is the 
(n + l)-bonacci constant. Hence, for k big enough, 

llfofelloo = Iflfel < aAia*-' 

The slowest growth, again modulo permutations, is obtained by those n- 
simplices corresponding to the bases {ei, e2+kei, • • • , e„-(-i-|-/cei}, whose barycen- 
ter bk = nkei + Y^^=i norm \\bk\\oc — nk + 1. □ 

Now we provide bounds for the volumes of the bodies Z{£i) in terms of the 
norms of the barycenters. 

Lemma 3. Let W — {wi,--- ,Wn+i) € M""*""^, n > 1, be a vector with non- 
negative components and let us build out of it a tree T{W) using the same 
algorithm used to build T{£), so that e.g. at the first tree level we find Wi = 
(wi, Wi + W2, ■ ■ ■ , wi + Wn+i) and the other n vectors obtained similarly. Then 
if the components of W satisfy the inequalities 

Wj <{n- l){wj^ + WjJ (1) 
the same inequalities hold for all other vectors of the tree. 

Proof. We prove the lemma by induction. Let us assume that the inequality 
is valid for all vectors up to the fc-th tree level and be W = (lu^, • • • ^w'^^j^i) 
one of the vectors at the level k. For the symmetry of the problem it is enough 
to verify that the inequality remains true for its first child W" — W[ and it is 
enough to check it in any two cases when its first component appears on the 
right side of the inequality and when it does not. 

In the first case let us assume ji = 1 and j2 = 2. Then the inequality reads 

n+1 

5^<<(n-l)K + u;g) 
i=3 

that is equivalent to (n — l)w[ + J2^=3 ^ {n — l)(2w[ + w'2) and therefore 

which holds by the inductive hypothesis. 
In the second case let us assume ji = 2 and j'2 = 3. Then the inequality 
reads 

«;i' + £<<(n-l)K+«;D 

J =4 
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that is equivalent to {'n~~l)w[+J2]=l ^'j — {n—l){2w[+w'2 + w';j) and therefore 

to w[ + X]j=4 5: ~ l)(^'i + ^'2 + ^3) + w'l which is true because, by the 
induction hypothesis, 

n+1 

w[ + Y^ w'j < (n - l){w'2 + W3) < (n - l){w[ + W2 + w'^) + w'l ■ 

All remaining inequalities are obtained by permuting the indices. □ 

Theorem 4. For every basis £ of K."+^ there exist real constants A, B such 
that 

A _ . B 



- < KZi) < 



for almost all multi-indices I . 

Proof. For this proof's sake it is convenient to use the same base £ and measure 
IJ, of Theorem [1] and the maximum norm for the barycenters. 

These choices have some important advantages: 1. if we call e^, j = 
1, . . . , n+1, the components of the vectors {ci} — £ with respect to the canonical 
basis of K"+\ then the vector W = (e"+\ . . . , e^^+J) built with the (n+l)-th co- 
ordinates of the basis vectors changes, when passing from the basis £1 to £iAk+it 
with the same rule illustrated in the Lemma above and satisfies the set of in- 
equalities H]); 2. if [/ii : • • • : hn+i] are the canonical homogeneous coordinates 
for MP", S{£) is entirely contained in the open set hn+i 7^ 0; 3. the component 
e""*"^ is not smaller than any other component for every i = 1, . . . ,n + 1; i. the 
expressions for the volume of S{£) and Z{£) are particularly simple. 

Now let {x], - ■ ■ , x"^^) be the components of the vectors of the basis £1, so 
that the homogeneous coordinates of the {n + 1) vertices of the n-simplex Sj 

"-I'll r,„A U^A^r 7. Ur, D . — T^l I ^1 . 



will be Ai — [xj : ■ ■ ■ : x"'^ ] and those of its body Zj will be B 



■ ^i'^^ + -^i^^]- ^ direct computation shows that 



1 



n+1 

n\ TT x"+i 



n 

i=l 



and 



where 6„ is the subdivision oi Zj in n + 1 simplices S'j*'', where each Sj'^^ has 
the same vertices of Si except for the fc-th vertex, which is replaced by the 
barycenter of Sj. 

Let us consider now one of the simplices s G 6„ and let [e^^ . ^ -I- 6^^. 3], j = 
1, . . . , n -|- 1, be its vertices. Note that, since all components of the basis vectors 
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are positive and no component is bigger than the last one, the barycenter's norm 
is 

n+l 

ii&/iioo-E^r+'- 
1=1 

Hence 



~ ' , - -^^ = 1 ^ ^ n+l , n+l ' - 



n+l 



Since we never used in our calculation the particular choice of the indices for 
the simplex s, these bounds are valid for all of them and therefore 

n + l + 



- < KZi) < 



n\\\bi\\^' -"^^ " - n\\\hi\\^^ 

□ 

Note that the inequality above does not hold for the n-simplices 5*/: for 
example, in the basis £ used above the simplices corresponding to the bases 

£k = {ei + fce„+i, • • • , e„ + fce„+i, e„+i} 

have barycenter 6fc = (1, • • • , 1, nfc + 1) and volume 

1 1 



• • • n\{k + l) 



which therefore is asymptotic to l/||6fc||" rather than to 

Numerical and analytical facts suggest that bodies' diameters are bound 
by the inverse of their barycenters' norm; in particular it is known to be true 
for n = 2 thanks to an indirect proof (see Section 14. 2[) and it is confirmed by 
numerical exploration of the rt = 3 case (see Section |473)) . We are led therefore 
to the following conjecture: 

Conjecture 1. For every basis £ ofM"~^^ there exist constants A and B such 
that 

A , , S 

< Z < 



\\b{Z)\\^ - ' ' - \\biZ)\\ 

where Z is any body associated to the tree T(£) and \Z\ its diameter with respect 
to the canonical distance d{[x\^ [y]) = arccos . 

As for the fractal dimension of F{£), we could not find any way to eval- 
uate exact non-trivial bounds for it; in next section we present the numerical 
evaluation of it for the cases n = 2, 3. 
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Figure 2: (a) log-log plot of the barycenters' norm vs. its index. The bi are arranged 
in the sequence naturally associated to the ordered tree T{£), namely bi follows bj if 
\I\ > \J\ or, in case the multi-indices have the same order, the lowest index which is 
different between I and J is bigger in I. Since there are 2* nodes at the level k the 
upper and lower bounds are evaluated using the fact that A: + 1 < ||i)2''-ill ^ ^^^2^^ 

and therefore log^k < \\bk\\ < ^(fc 1)'°^= (b) log-log plot of the length of the 
1-simplices Si vs. the barycenters' norm. 

4 Analysis of the cases n = 1, 2, 3 
4.1 The case n=l 

The construction we discussed above does not strictly speaking apply to the 
n = 1 case. E.g. bodies here are simply single points, Theorem [3] does not 
apply and all asymptotics about the measures of bodies have no meaning here. 
Nevertheless a iew things survive: the tree T{£) and its Fibonacci sections can 
still be built and we can study the asymptotics of the lengths of the 1-simplices. 
To begin, let us choose 

£ = {ei = (l,0),e2 = (l,l)} 

The set F{£) is invariant with respect to the projective transformations 

M[hi ■■ h2]) = [hi + h2 : hi],M[hi ■ h2]) = [2/ii - /12 : hi] 

and it is obtained from the segment [0, 1] (in the projective chart x = 1) by 
removing a countable set of infinite (rational) points, so that it has full measure 
and therefore dim/f F{£) = 1. 

The growth rate of the Fibonacci sections here is given by the Golden Ratio 
a = (1 -I- V5)/2 and for the norm of the sections' barycenters we have the 
inequalities k + 2 < \\bk\\oo for the slowest section and ||&fe|joo < ^'^'^^^ ^'^^ 
fastest (see fig. [21[a)). In particular the components of the two root Fibonacci 
sections are exactly the Fibonacci numbers: e.g. taking bi = h{Z£) = (1,2) and 
62 = 6(2'{ei, 61-1-62}) = (2,3) we have that 63 = (3,5), 64 = (5,8) and so on. 

Asymptotics of bodies have no meaning here but still we can say something 
about the asymptotics of the lengths of the 1-simplices constituting the binary 
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tree T{£). Indeed if £' = {aei + be2, cei ^ 
chart y — I, 

fji[E') = d{[aei + 662], [cei + 1^62]) 
\ad - bc\ ■ \x1y2 - X2yi\ _ 



de2} with Ci = {xi,yi) then, in the 
axi + bx2 cxi + dx2 



1 



ayi + by2 cyi 



dy2 
1 



(ayi + by2){cyi + dy2) {ayi + by2){cyi + dy2) {a + b){c + d) 

where \x1y2 — X2yi \ = 1 is the surface of the parallelogram corresponding to £ 
and \ad — bc\ = 1 because of the way the algorithm produces the new bases. In 
our concrete case yi — 1 and therefore 



-(- 



1 



1 



< 



\\b\\io ' \\b\\oo'a + b • c + d' - \\b\\ 

Numerical illustrations of this pair of inequalities are shown in fig[2l^b). 



4.2 The case n=2 

This is the only case where the polytopes corresponding to the bases and to the 
bodies are of the same kind, namely triangles. The algorithm that produces the 
fractal reduces here to the following: 

Algorithm 1. 

1. On the three edges of the triangle A with vertices {[ei]}i=i_2.3 select the three 
points fi = [62 + 63], /2 = [es + ei], /s = [ei + 62]; 

2. subtract from A the interior of the triangle Z (the "body" of A) with vertices 
{fli f2, /a};' 

3. repeat recursively the algorithm on each of the three triangles that are left 
after the subtraction. 

Note that no two bodies have in common more than a point, i.e. they 
meet transversally, so the set F{£) is never empty and actually it contains 
uncountably many points; countably many of them can be explicitly evaluated 
through Fibonacci sections of the ternary tree T{£). 

Consider for example the case 

£^{ei = (1, 0, 0), 62 = (0, 1, 0), 63 = (0, 0, 1)} . 

As shown in Theorem [3] the barycenter of every body triangle is the vector 
sum of the barycenters of the three body triangles it touches with its vertices 
(note that by construction no two bodies have a vertex in common) and when a 
body touches one of the sides of the root triangle S{£) then we sum instead the 
vectors (—1,1,1), (1,-1,1) and (1,1,-1) in correspondence respectively with 
the sides opposed to the vertices [ci], [62] and [63]. The barycenters of one of 
the six root Fibonacci sections are determined by the first elements 



6_3 = (1,1,-1), ^^-2 = (1,-1,1), 6-1 = (-1,1,1) 
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so that the generic element of the section is given by bk = (a/c_2, a^-i, Ofc), 
where Uk is the sequence of Tribonacci numbers with initial conditions a_2 = 1, 
a_i — 1, flo = 1. The expression of the generic term is given by 

_ k I (l-a)(l-/3) 

where a, p and /3 are the roots of the Tribonacci equation ~ x'^ + x + 1. 
Since |/3| < a we have that ||6fe||oo < 3^2*^2^-1 '-'''' ^^"^ hmit point (see fig. [3]) 
is (1 : a : a^). Note that all barycenters of this sequence lie on the "projective 
Tribonacci spiral" 

j{t) = [a{t-2) : a{t~ 1) : a{t)] 

where a{t) is the trivial analytical extension of the Uk sequence. The fractal is 
invariant with respect to the projective transformations 

^i{[hi : /i2 : hs]) =[hi + /i2 + : /12 : 
V'2([/ii : h2 ■■ hs]) =[hi : hi + h2 + h3 : /13] 
'ipsiihi : /i2 : ^3]) =[^1 : h2 ■■ hi + h2 + h^] 

so by applying any finite composition of them we obtain countably many explicit 
points of F{£). 

The slowest sections in the barycenters' norms growth is, modulo indices 
permutations, 

ifc = . . . ^ 1 = {e^'"'' = ei, 62''^ = 62 + fcei, e^*"^ = 63 + kei} 

k 

for which hk — {2k + 3, 1, 1) and therefore ||&fc||oo — 2k + 3. 

In figs. I4.2r a-c) we show the numerical results for the asymptotic behaviour 
of the barycentric norms and the bodies' surfaces and diameters. 

Note that in this particular case Conjecture[T]is known to be true through an 
indirect proof. Indeed, this fractal comes up naturally in the study of the asymp- 
totics of plane sections of periodic surfaces, which in turn comes from the prob- 
lem of the motion of quasi-electrons under a strong magnetic field (see |NM03j 
for a detailed account), in the particular case of the regular triply-periodic skew 
polyhedron {4,6|4} |DD09| . In that setting the basis is 

£c = {ei - (1, 0, 1), 62 = (0, 1, 1), 63 = (1, 1, 0)} 

and the barycenter 5 of a body Z represents a homological discrete "first inte- 
gral" of a Poisson dynamical system which dictates the asymptotic directions 
of the plane sections in the following way: the open sections obtained by cut- 
ting the polyhedron with planes perpendicular to every direction a; G Z are all 
strongly asymptotic to the direction "w x 6". It is a general theorem of that 
theory the fact that the diameter of a body Z is bounded by C/|j5(Z)|| where 
C is a constant depending only on the surface [De 05j . which then establishes 
the following theorem for this n = 2 case: 
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(0:1:1) 




(1:0:3) (1 . : 2)(2^o^3) (1:0:1) 



Figure 3: Plot of F^{£), namely of the bodies up to the forth recursion level, for 
£ = {(1,0, 1), (0, 1, 1), (0,0, 1)} in the = 1 projective chart of RP^. Bodies are 
colored in green, so the points of are the white ones. The homogeneous coordi- 
nates of the vertices of the first and second level bodies are shown together with the 
body's barycenters, for which we used the square brakets for sake of clarity. Note that 
barycenters can be obtained in three ways: (i) by summing the barycenters of the three 
bodies touched by the vertices - note that in case a vertex touches a root simplex edge 
then the following should be used: (1, 1, 1) for the edge opposite to [0 : : 1], (1, —1, 1) 
for the one opposite to [1 : : 1] and (—1, 1, 1) for the one opposite to [0 : 1 : 1]; (ii) 
by summing the coordinates of the vertices of the triangle that generated the body; 
(iii) by summing the coordinates of the vertices of the body and dividing them by 2 
- this corresponds to the fact the volume associated to the basis corresponding to the 
vectors {ei + 62, 62 + 63, 63 + ei} is double with respect to the basis {ei, 62, 63}. 
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Figure 4: Plot of the barycenters of all 3487590 bodies up to the thirteenth recursion 
level for £ = {(1,0, 1), (0, 1, 1), (0,0, 1)} in the = 1 projective chart of RP^. The 
colors of the points goes from red to blue as the Euclidean norm of the barycenters 
grows. By Corollary [1] the closure of the set of barycenters is equal to the bound- 
aries of all bodies plus the points of the fractal F{£), so this picture represents an 
approximation of the real fractal (and actually no point shown belongs to F{£) since 
barycenters are all contained inside the bodies. Nevertheless they can approximate as 
close as wished the set F{£) and so they can be used to derive a numerical evaluation 
of the box-counting dimension of F{£). 
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Theorem 5. Let Z be a body in T{£c) with area ii{Z) (where fi is the same 
measure used in Theorem^, diameter \Z\ and bary center b. Then the following 
inequalities hold asymptotically: 

1 1^1 6 1 12\/3 

Proof. The inequality for the area of Z is just the restriction of Theorem [4] to 
n = 2 together with the fact that ||6||oo < \\b\\ < \/3||&||oo- The right hand 
side for the diameter comes from the general theory of plane sections of a triply 
periodic surface that, apphed to this particular case, states [De 05j that the 
distance between the barycenter and the bodies' vertices is bounded by 3/||6||, 
where the 3 is the double of the area of the basic cell of the periodic surface 
cited above in this section. The left hand side comes simply from the fact that 

a triangle of area a cannot have a diameter smaller than '\J 2a/\/3. □ 

Being unable to evaluate analytical bounds for the Hausdorff dimension dc of 
F{£c), we compute numerically four different quantities that may give hints on 
whether dc is integer or not (the non-integrality of dc would confirm a general 
conjecture by Novikov [NMOSj ). 

First we get a direct upper bound for the Hausdorff dimension by counting 
the smallest number of squares of side e = 2~', / = 0, • • • , 12, needed to cover 
F^^, i.e. the union of all bodies up to the 12-th order of recursion; as shown in 
fig.[12i;h), we get dc < 1.7. 

Then we evaluate the Minkowsky dimension, namely the limit 

2_,^logm) 
£-^0 log e 

where V{F^) is the surface of the e neighborhood of F, using the formula [Fal97[ 
IGaiOBj 

Ale Ale A;e 2 

V{F,)^pe + eJ2p^ + A~Y.''^ + ''(^-J2i^) 

i—l 2—1 i—1 

where is the integer such that Pk,+i ^ e ^ Pk^^ Pk is the radius of the inscribed 
circle to the body Zk and the bodies are sorted in descending order with respect 
to the radii. In fig. I4.2r g) we show the numerical results we got by evaluating 
the volume of the neighborhoods of £ of radii r„ — 1.2~" for n = 1, • • • , 50, 
which suggests a Minkowsky dimension between 1.7 and 1.8. 

Next, we evaluate numerically the growth rate of the radii after sorting 
them in decreasing order (fig. I4.2f d)) and then the corresponding bounds for 
the bodies areas ffig. I4.2f e)) and diameters ffig. I4.2f f)). In this case we obtain 
that e X fc--69, Ak-'^-^^ < ak < Bk-^-^ and A'k' ''^ < Pk < B'k" '^, so that 
< < iJ"e-i45^ From this we get a second evaluation, compatible but 
much looser, for the Minkowsky dimension: 1.35 < dimM^"(i^^c) < 1.86. 

Finally, we use Corollary [1] and evaluate the box-counting dimension of the 
set B ~ {[bi{£c)]} C RP^, namely the set of barycenters of the bases in the 
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Figure 5: Log- log plots for the main quantities in tlie n = 2 case for £ = 
{(1, 0, 1), (0, 1, 1), (0, 0, 1)} (a) Barycenters norms vs indices - as explained in fig. [5] 
the bk are ordered according to the natural order induced by the tree, so that 
2fc + 3 < 1163^ II < V3aili5|liizMat = Aal and therefore log fc + logg 12 < 

||6fe|| < j4fc'°^3"^. (b) Bodies' volumes vs barycenters norms and (c) bodies' diameters 
vs barycenters norms - the lines bounding the numerical data come immediately from 
the inequalities in Theorem [5] For the next three plots no exact formulae are known 
so the lines shown represent just an interpolation of the numerical data, (d) Radii of 
the circles inscribed in the bodies vs k after sorting the radii in descending order, (e) 
Areas of the bodies and (f) their diameter sorted according with their radii. 
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tree T{£c) considered as points in the projective plane. Since the closure of 
B is the union of F{£c) with the (one-dimensional) boundaries of the bodies, 
a dimension higher than one must be due to the points in F{£c)- As shown 
in Fig. 14. 3| the dimension appear to be about 1.69. Notice that applying this 
method to the Sierpinki triangle, whose dimension is ds = log 3/ log 2, gives the 
correct approximation to the third digit ds — 1.59. 

In conclusion, the four evaluations are in excellent agreement with each other 
and indicate a non-integer Hausdorff dimension for this fractal, probably about 
1.7; finding exact bounds would be nicer though since this would represent the 
first analytical confirmation of a conjecture of Novikov about the non-integer 
dimension of fractals coming from the theory of asymptotics of plane sections 
of triply-periodic surfaces. 



4.3 The case n=3 

When n = 3, every body has 6 vertices: one for each edge of the tetradedron 
they belong to, and eight triangular faces, one for each face and one for each 
corner of the tetrahedron. Bodies that touch each other share a whole triangle 
(rather than a single point as in the n = 2 case) in the following way: bodies 
can meet only on the faces that do not come from the tetrahedra Si and, on 
those faces, these shared triangles form a fractal of the n — 2 kind (see Fig. l4.3p . 

In the particular case of £t = {(1, 0, 0, 1), (0, 1, 0, 1), (0, 0, 1, 1), (0, 0, 0, 1)} 
the barycenter of the root tetrahedron is (1,1,1,4) and the volume inequalities 
translate in 

2 2 • S'^ 

< /i(Z/(fT)) < 



4 



In case of barycenters' norms we have 4 + 3fc < for the slowest tree section 
and llfoJI < 2^-^^§Mii^lLQ,fe+4 ^j^g fastest. 

II '^11 — (a4-/3)(a4-/3)(a4-7) ^ 

Numerical evaluations of the Hausdorff dimension are more cumbersome for 
n = 3 because the number of bodies grows very large after few iterations of the 
generating algorithm (getting rather heavy on both CPU and RAM consump- 
tion) and their geometry gets much more complicated. 

First we evaluate the Minkowsky dimension as the growth rate of the volume 
V and surface 5* of the bodies when sorted by the radius p = V/S. In this case 
we use the fact that |Gai06| 

i=k, + l i=l i=fcj + l 

where is the volume of the neighborhood of F{£t) of radius e, the integer 
such that Pk,+i < e < Pfe, , Vi and Si the volume and surface of the body Z;, 
S and H the surface and mean curvature of the starting tetrahedron. From 
the numerical data (see fig. I4.3r a-c)) we obtain that e x fc~-^^, Ak^-'''^ < Vk < 
Bk'^-^ and < Sk < B'k--^, so that A"e^-^^ < < B"e-^^^ and 

therefore 1.66 < diuiM F < 2.75. Unfortunately, milikc in the n = 2 case, we 
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are not able to exclude from this bounds that the fractal has integer dimension 
dimM F = 2; note that this is exactly what happens for the Tetrix, i.e. the 
three-dimensional analog of the Sierpinki triangle. 

Next we use Corollary [T] and evaluate numericlaly the box-counting dimen- 
sion dbc of the set B = C MP"^, namely the set of barycenters of the 
tree T{£t) considered as points in the projective three-space. Analogously to the 
case n — 2, the closure of B is the union of F{£t) with the (two-dimensional) 
boundaries of the bodies, a dimension higher than two must be due to the 
points in F{£t)- We obtain d^c — 2.20 (see Fig. 14. 3|) . Notice that the very same 
method, applied to the Tetrix, whose dimension is — 2, gives the quite close 
result dr ^ 2.01. 

We could not get useful information from the other two methods used in the 
n — 2 case. In conclusion, the two numerical results we obtained are compatible 
with each other and the estimate of the box-counting dimension of barycenters 
is sufficiently far from integer to make us think that, unlike the Tetrix, this 
fractal may have Hausdorff dimension higher than 2. 
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Figure 6: log- log plots for the main quantities in the n — 3 case for £t ~ 
{(1, 0, 0, 1), (0, 1, 0, 1), (0, 0, 1, 1), (0, 0, 0, 1)} (a) Barycenters norms vs indices - as ex- 
plained in fig. [2] the bt are ordered according to the natural order induced by the 
tree, so that 3k + i < ||6,fe ,11 < 2a|-^i^4#-^|4ii^^Qt = Aa'l and therefore 

log 4- log4 108 < ||6fc|| < Ak^°^'^°"^. (b) Bodies' volumes vs barycenters norms - 
the lines bounding the numerical data come immediately from the inequalities in The- 
orem [5] (c) Bodies' diameters vs barycenters norms and (d) bodies' surfaces areas; 
for these quantities no exact formulae are known so the lines shown represent just an 
interpolation of the numerical data. 
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Figure 7: Comparison of tlie numerical data between tlie n — 3 case F{£t) for 
£t = {(1,0,0, 1), (0,1, 0,1), (0,0, 1,1), (0,0, 0,1)} and the Sierpinski tetraliedron S 
about the asymptotic behaviour of the bodies' volumes V and surfaces S sorted by 
their "radii" p — V/S in descending order. No exact formulae are known for these plots 
so the lines shown above represent just an interpolation of the numerical data. (a,c) 
Radii of the circles inscribed in the bodies vs k after sorting the radii in descending 
order for the n = 3 case (left) and the Sierpinski tetrahedron (right). (b,d) Areas of 
the bodies and (c,f) their surface sorted according with their radii for the n = 3 case 
(left) and the Sierpinski tetrahedron (right). From these interpolation we obtain that 
for F{£t) we have 1.6 < dims F{£t) < 2.9 and for the Sierpinski tetrahedron we get 
the correct answer dim 5 = 2. 
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Figure 8: Plots relative to the evaluation of fractal dimensions of Fc = F{£c) C RP^ 
(see Section 1121 and Ft = F{£t) C RP^ (see Section g^J. (a) Evaluation of 
the Minkowsky dimension by direct numerical computation of the area of the of 
the e neighborhood of Fc, based on the fact that, for "nice" fractals F, dim_F = 

limj_,o+ 1^2— '"og^' j • (b) Evaluation of the box-counting dimension dbdFc) by di- 
rect computation of the number of squares needed to cover the 12-th order 
approximation of Fc, for e = 2"*^, k = 1, • • • ,12. The data strongly suggest that 
dbc{Fc) — 1.7. (c) Evaluation of the box-counting dimension d'fj^Fc) of the set of 
barycenters (considered as points in RP'^) of the bodies Z{£j), £j £ T{£c)- By Corol- 
lary [TJ the closure of this set is the union of F{£c) with a 1-dim. set. We get, in 
excellent agreement with (b), ci'(,c(^c) — 1-69. (d) Same evaluation as in (c) in case 
of the Sierpinki triangle. The estimated dimension dj,^, ~ 1.59 is in perfect agreement 
with the exact result log 3/log2 ~ 1.585. (e) Evaluation of the box-counting dimension 
d-bci^T) of the set of barycenters (considered as points in RP^) of the bodies Z{£i), 
£i £ T{£t)- By Corollary [TJ the closure of this set is the union of F{£t) with a 
2-dim. set. We get d'^dFT) — 2.20. (f) Same evaluation as in (e) in case of the Tetrix 
(three-dim. analogue of the Sierpinki triangle). The estimated dimension d'^^ ~ 2.01 
is in excellent agreement with the exact result 2. 
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(c) (d) 

Figure 9: (a,c) Total view and detail of F^{£), namely of the bodies up to the 
fifth recursion level, for £ = {(1, 0, 0, 1), (0, 1, 0, 1), (0, 0, 1, 1), (0, 0, 0, 1)} in the /i" = 1 
projective chart of RP^. The bodies of S{£) are shown, up to the fifth recursion level, 
in red-green colors; the points of F^{£) are their complement in the tetrahedron of 
vertices (in the chart /i* = 0) (1,0,0), (0,1,0) and (0,0,1). (b,d) Total view and 
detail, up to the forth recursion level, of the Sierpinski tetrahedron. 
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